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Understanding the acceleration of the universe and its cause is one of the key problems in physics and cosmology 
today, and is best studied using a variety of mutually complementary approaches. Daly and Djorgovski (2003, 
2004) proposed a model independent approach to determine the expansion and acceleration history of the universe 
and a number of important physical parameters of the dark energy as functions of redshift directly from the data. 
Here, we apply the method to explicitly determine the first and second derivatives of the coordinate distance with 
respect to redshift, y' and y" , and combine them to solve for the kinetic and potential energy density of the dark 
energy as functions of redshift, K{z) and V{z). 

A data set of 228 supernova and 20 radio galaxy measurements with redshifts from zero to 1.79 is used for this 
study. Values of y' and y" are combined to study the dimensionless acceleration rate of the universe as a function 
of redshift, q[z). The only assumptions underlying our determination of q(z) are that the universe is described 
by a Robertson- Walker (RW) metric and is spatially flat. We find that the universe is accelerating today, and 
was decelerating in the recent past. The transition from acceleration to deceleration occurs at a redshift of about 
zt = 0.42±o;o|. Values of y' and y" are combined to determine K{z) and V{z). These are shown to be consistent 
with the values expected in a standard Lambda Cold Dark Matter (LCDM) model. 



The acceleration of the universe at the present 
epoch has been studied in the contexts of specific 
models using coordinate distances to type la su- 
pernovae |1I2I3I4I5I6I7| . and FRII radio galaxies 
|8l9ll0lll] . in addition to other techniques. These 
studies indicate that the universe is expanding at 
an accelerting rate at the present epoch. Gen- 
erally, these studies are done in the context of a 
specific cosmological model, such as an open uni- 
verse with non-relativistic matter, a cosmological 
constant, and space curvature (e.g. T"2"8''9'), a 
spatially flat universe with non-relativistic matter 
and dark energy that has an energy density that 
can evolve with redshift but which maintains a 
constant equation of state (e.g. 2 3 4 5 6.7.10, '). 
or a spatially flat universe with non-relativistic 
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matter and an evolving scalar field (e.g. |12lll| '). 

In each of these studies it is implicitly assumed 
that the universe is desribed by a RW metric and 
that General Relativity is the correct theory of 
gravity; in addition, a particular functional form 
for the redshift evolution of the energy density 
of some new component is assumed. The data 
are then used to constrain the parameters that 
describe the assumed functional form for the red- 
shift evolution of whatever was being considered 
as the driver of the acceleration of the universe. 

A complementary approach was suggested by 
|13ll4j who showed that the recent expansion and 
acceleration history of the universe, and some 
properties of the driver of the accelertion, can be 
determined directly from the data after specify- 
ing a minimal number of assumptions. Assuming 
only that the universe is described by a RW met- 
ric and is spatially flat, the data can be used to 
solve for the dimensionless expansion and accel- 
eration rates of the universe as functions of red- 
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Figure 1. Dimensionless coordinate distance y(z) 
to 71 Legacy and 157 Gold supernovae, and 20 
radio galaxies; y{z) = Ho{aQr), Hq is Hubble's 
constant, and (aor) is the coordinate distance to 
a source at redshift z. It is convenient to work 
with y{z) because it is independent of Hq (i.e. 
(aor) oc Hq^, so HQ{aor) is indepenent of Hq). 
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Figure 2. Results obtained with the mock data 
set of 248 sources described in the text. The re- 
sults are in an excellent agreement with the input 
cosmolgoy, with no apparent bias. 
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Figure 3. As in Fig. |21 for y". The correct 
assumed cosmology is recovered with a negligible 
bias. 
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Figure 4. The first derivative of the coordinate 
distance with respect to redshift for the actual 
data set of 248 sources. The zero redshift value we 
measure is y^ = 1.025±0.022; the predicted value 
in all models is 1.000. The values for the stan- 
dard LCDM model with = 0.7 and flom = 0.3 
are shown as the solid line in this and all sub- 
sequent plots. Best fit Cardassian (dotted line) 
and Chaplygin Gas (dashed line) models are also 
shown, and are described in the text. 
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Figure 5. The second derivative of the coordi- 
nate distance with respect to redshift for the ac- 
tual data set of 248 sources. The measured zero 
redshift value is i/'q = —0.55 ±0.10; the value pre- 
dicted in the LCDM model shown is 0.45. 
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Figure 6. The deceleration parameter qiz), 
where q{z) = -[1 + + z)ly'\ 13 . The zero 
redshift value is go = -0.46 ± 0.08. The pre- 
dicted value in the LCDM model shown is —0.55. 
Our fits are systematically higher than the LCDM 
model shown by about Icr. 



shift, E(z) and q(z), respectively. This can be 
done without specifying a theory of gravity, or 
anything else. The function q{z) thus obtained is 
a direct measure of the acceleration/deceleration 
of the universe at different epochs. The key ingre- 
dients that go into the determination of E{z) and 
q[z) are the first and second derivatives of the co- 
ordinate distance with respect to redshift, dy/dz 
(or y') and d'^y/dz^ (or y"), which are obtained 
from the coordinate distances to supernovae and 
radio galaxies at known redshift, as described by 
|13ll4j . Thus, rather than assuming a functional 
form for the redshift evolution of the "dark en- 
ergy" and constraining the model parameters, it 
is possible to solve for quantities such as q{z) di- 
rectly. 

This direct approach indicates that the uni- 
verse is accelerating today, and was decelerating 
in the recent past. The data used for the results 
shown here include the sample of 157 "Gold" su- 
pernovae 0, the sample of 71 supernvae from the 
Supernova Legacy Survey U], and the 20 radio 
galaxies of IH], as described in detail ^Hl- The 
total sample of 248 sources is shown in Fig. ^ 
there are no systematic differences seen among 
the three groups of measurements in the redshift 
ranges of their overlaps. 

The first and second derivatives of the coor- 
dinate distance with respect to redshift are ob- 
tained using the numerical differentiation method 
described by |13I14| . To test whether the method 
introduces a bias in the results, a mock data set 
of 248 sources with the same redshift distribu- 
tion and fractional uncertainty per point as the 
actual data was constructed assuming a LCDM 
model with fiom = 0.3 and J^a = 0.7, and ana- 
lyzed. The results are shown in Figs. |2 and [3]. 
We see that no bias has been introduced by the 
numerical differentiation technique. 

Values of y' and y" are shown on Figs. @] and 
|S] . The ringing seen in these figures is most likely 
due to sparse sampling. In these plots, and in the 
ones that follow, we do not consider these fluctua- 
tions at higher redshifts to be statistically signifi- 
cant, as they are commensurate with our derived 
l-cr error bars. The results are consistent with 
the LCDM model. The LCDM model is based 
on General Relativity with non-relativistic mat- 
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Figure 7. The kinetic energy of the dark energy 
K(z) in units of the critical density today. The 
zero redshift value is Kg = 0.02 ± 0.03; the ex- 
pected value for the LCDM model shown 0. 
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Figure 8. The potential energy of the dark en- 
ergy V{z) in units of the critical density today. 
The zero redshift value is Vq = 0.62 ± 0.05; the 
expected value for the LCDM model shown is 0.7. 



ter f^om and a cosmological constant. This pro- 
vides an excellent description of the data. Curves 
showing predictions in two modified gravity mod- 
els in a spatially flat universe are also shown on 
Fig. These curves shown are obtained using 
the best fit model parameters obtained by ^Hj for 
the Cardassian model of and the generalized 
Chaplygin gas model of JH| based on the model of 
[TO] ; this is consistent with the results obtained by 
|20| . Clearly, the LCDM model provides a better 
description of the data than do either of the mod- 
ified gravity models. Thus, this large-scale test of 
General Relativity shows that GR provides an ex- 
cellent description of the data on these very large 
length scales of about 10 billion light years. 

The deceleration parameter q(z) is shown on 
Fig- El- These results allow a determination of 
the redshift at which the universe transitions from 
an accelerating phase to a decelerating phase; we 
find this redshift to be zt = 0.42±g Q|, consistent 
with the values quoted by |B] and |13I14| . The up- 
per bound on this transition redshift is uncertain 
because of the fluctuations in q{z) which arc due 
to sparse sampling at high redshift. 

It is well known that K = 0.5{p + P) and 
V = 0.5(p — P), where p and P are the energy 
density and pressure of the dark energy. In |14| 
we show that both p and P may be written in 
terms of the first and second derivatives of the co- 
ordinate distance. Combining these, we find that 
{K/poc) = -(1 + z)(y")(y')"V3-0.5^!nm(l+^)^ 
and (V/poc) - {y')-^[l + il + z)y"{y')-'m - 
0.5Slom(l + z)^ , where poc is the critical den- 
sity at the current epoch. These are shown in 
Figs. |Zj and |H]. In obtaining K and V, the as- 
sumptions made to obtain P and p apply: the 
universe is spatially flat; the kinematics of the 
universe are accurately described by general rela- 
tivity; and two components, the dark energy and 
non-relativisitc matter (with flom = 0.3), are suf- 
ficient to account for the kinematics of the uni- 
verse out to redshift of about 2 (see the discussion 
in Functional forms for P(z) and p{z) for 

the dark energy are not assumed, nor is any as- 
sumtion made regarding the equation of state of 
the dark energy. The work presented here on the 
potential energy, V{z), is complementary to the 
work of^mMMM- 
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Thus, our (nearly) model-independent method 
provides results which are consistent with those 
from the more traditional approaches, in a largely 
complementary fashion; at the very least, it is a 
new way of looking at the data. As the quality 
and size of relevant data sets increase, we can ex- 
pect even more useful constraints to emerge from 
this approach. 
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